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Abstract 

We propose a scheme to assign a 4-momentum to each WIMP in new physics event producing 
a pair of mother particles each of which decays to an invisible weakly interacting massive particle 
(WIMP) plus some visible particle(s). The transverse components are given by the value that 
determines the event variable Mt2 , while the longitudinal component is determined by the on-shell 
condition on the mother particle. Although it does not give the true WIMP momentum in general, 
this Mr2-assisted on-shell reconstruction of missing momenta provides kinematic variables well 
correlated to the true WIMP momentum, and thus can be useful for an experimental determination 
of new particle properties. We apply this scheme to some processes to measure the mother particle 
spin, and find that spin determination is possible even without a good knowledge of the new particle 
masses. 
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1 



The Large Hadron Collider (LHC) at CERN will explore soon the TeV energy scale where 
new physics beyond the Standard Model (SM) is likely to reveal itself. There are two major 
motivations for new physics at the TeV scale, one is the hierarchy problem and the other 
is the existence of dark matter (DM). Constraints from electroweak precision measurements 
and proton decay suggest that TeV scale new physics preserves a parity under which 



the new particles are odd, while the SM particles are even. We 
the weak scale supersymmetry (SUSY) with conserved i?-parity 



1-known examples include 
1|, little Higgs model with 



T-parity Q|, and extra-dimensional model with KfC-parity j^. The lightest new particle in 
these models is typically a weakly interacting massive particle (WIMP) which is a good DM 
candidate. 

A typical LHC signature of new physics model with conserved parity is the production 
of a pair of new particles decaying to invisible WIMPs plus some visible SM particles. As 
the WIMP momenta can not be measured, an event by event reconstruction of the new 
particle's rest frame is not available, which makes the determination of new particle mass 
and spin quite non-trivial. For the mass measurement, several methods have been proposed 
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so far [J, 15|, 16[ , and they might work with certain accuracy depending upon the kinematic 
situation. One can also attempt to determine the spin from the production rate or from 
a kinematic variable distribution that shows spin correlation However, usually spin 
determination is more difficult as it requires a larger statistics and/or a good knowledge of 
the mass spectrum and branching ratios. 

In this paper, we propose a scheme to assign a 4-momentum to each WIMP in new physics 
event, which can provide kinematic variables useful for the experimental determination of 
new particle properties, particularly the spin. We dub this scheme the " MT2-Assisted On- 
Shell" (MAOS) reconstruction as the transverse components are given by the value that 
determines the collider variable Mj-i [sl, while the longitudinal component is determined by 
the on-shell condition imposed on the mother particle. This MAOS reconstruction of WIMP 
momenta can be done even without a good knowledge of the WIMP and mother particle 
masses. In the following, we discuss the MAOS reconstruction in more detail, and apply it 
to the 3-body decay of gluino and the 2-body decay of Drell-Yan produced slepton, as well 
as to their universal extra-dimension (UED) equivalents, in order to see the spin effects in 
these processes js, ol. 
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To start with, let us consider the new physics process: 

pp + F(2) ^ v{p,)x{h) + V{p2)x{k2), (1) 

where V{p) denotes generic set of visible SM particles with total 4-momentum p^ and x(^) 
is the WIMP with 4-momentum fc^. For an event set of this type, one can introduce trial 
mother particle and WIMP masses, my and m^, and impose the on-shell condition together 
with the missing Et constraint: 

{p, + k,Y = ml, k1 = m% klT + k2T = P^^ (2) 

where ^''t^^^ denotes the missing transverse momentum of the event. As these provide 6 
constraints for 8 unknowns, /cf {i = 1,2), there are 2-parameter family of solutions, which 
can be parameterized by ki^. For any choice of real ki^, e.g. ki^ = k^, the longitudinal 
WIMP momenta are determined to be 
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^ Kl, (3) 



where E^t = \/pI + |PiTp, E}j, = y^m^ + Ik^rP, and Ai = \{m\ -m^,- pf) + p^^ • k^j. for 
kiT = kr and k2T = p™*^*^ — kr- It is obvious that kf^^ are real if and only if \Ai\ > E^Efj, 
which is equivalent to 

my > max{MP , M^'^ } , (4) 



where = Jpf + + 2( E^j^Efj, — piT ■ kjr) corresponds to the transverse mass of the 



mother particle Y{i) with kij- = k^- and k2T = p™'^^ — kj^. 

In principle, one could choose event-by-event any value of k^^. However, the condition (jl]) 
suggests that the best choice of ky is the one minimizing ma:x{M^\M^^} for each event, 
i.e. the value giving the collider variable Mt2 SJ: 

MT2{Pi,m^)= min [max{M^^\ M^^^}] , (5) 



where the minimization is performed over all possible WIMP transverse momenta kj^ under 
the constraint ki^ + k2T = p™^*^*^. For given trial masses ^^^^ choice of k^, which is 

unique for each event, allows the largest event set to have real kfj^. In the following, we call 
this scheme the Mr2-Assisted On-Shell (MAOS) reconstruction, which assigns one or both 
of the two 4-momenta 



fc± = (Jm2 + |kT|2 + |fc±|2,kT,A;±), (6) 



to each WIMP in the new physics event ([T])^. With these MAOS momenta, one can con- 
struct various kinematic variables whose distribution shape may provide information on new 
particle properties. 

By construction, if my > M^^^lm^), where M^2^ denotes the maximum of Mt2 over 
the event set under consideration, the reconstructed MAOS momenta are real for all events. 
However, if my < M^2^{m~^), generically there are some events which do not allow real 
longitudinal MAOS momenta. Also by construction, for the Mt2 end-point events with 
= M^2^, the MAOS momentum is same as the true WIMP momentum if (i) the trial 
masses used for reconstruction are same as the true masses, and (ii) the considered event set 
is large enough to give M™2^''(m^ = m^™^) = my'^'^. (Note that M™2^''(m^ = m^^'^^) < my^^ 
for generic event (sub) set, and the bound is saturated if the full event set is taken into 
account.) One can also find that = for an event whose Mr2- value is same as my, 
which obviously includes the end-point event when the trial mother particle mass is chosen 
as my = M™2^{m^. In fact, the most interesting feature of the MAOS momenta, which 
will be discussed below with a specific example, is that they are distributed around the true 
WIMP momentum with non-trivial correlation even when constructed with generic trial 
WIMP and mother particle masses. 

As a specific application of the MAOS reconstruction, we first consider the symmetric 
3-body decays of gluino pair in SUSY model, g + g ^ qqx + QQXi ^ind also the similar 
decays of KK gluon pair in UED-like model. For SUSY case, we choose the focus (SPS2) 
point of mSUGRA scenario which gives the following weak scale masses: m|''™ = 779 GeV, 
^true _ ]^22 GeV, and mj^™*' ~ 1.5 TeV. We also consider its UED equivalent in which 
the gluino is replaced with the first KK gluon (yf(i), the Bino B with the first KK U{1)y 
boson B{i) , and squarks with the first KK quarks. Using MadGraph/MadEvent 1^, we have 
generated the events at parton-level for both SUSY and UED cases, and constructed the 
MAOS momentum /c^. We then examined the distribution of AA; = A;^ — for both SUSY 
and UED event sets with various choices of trial masses m^ y and also of the event subset 
selected by their Mt2 values. We found that the distribution of Akx is always peaked at zero 
for generic value of m^, and its width gets narrower if one considers an event subset including 



^ For simplicity, here we consider the production of YY in which the up-streaming momentum is neghgible. 
It is however straightforward to generahze the MAOS reconstruction to the process with a sizable up- 
streaming momentum. 
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only the near end-point events of Mt2- On the other hand, the distribution of Afc/, is more 
chaotic, partly because of the error propagation from Akx and also the two-fold degeneracy 
of the longitudinal component. Still it is peaked at zero, although the width is significantly 
broader, for a wide range of (m^,my) which includes the case with my = M^2^{m^. As 
an example, we depict in Fig. [1] the distributions of AkT,L (including both of kf^) for the 
MAOS momenta of the SPS2 SUSY event set, which has been constructed with = and 
my = M^2^{m^ = 0). Fig. [TK shows the distributions of the full event set, while Fig. [lb is 
for a subset including only the top 10% end-point events of Mt2- We can see that the MAOS 
momentum has a good correlation with the true momentum even for the full event set, and 
the correlation becomes stronger for the near end-point events of Mt2- This suggests that if 
one has an enough statistics, it can be more efficient to do MAOS reconstruction using only 
the near end-point events. 




FIG. 1: The distributions of — k^"^^^ for (a) the full event set, and (b) the top 10% end-point 
events of Mt2- Here the MAOS momenta were constructed with = and my = M^2^{m^ = 0). 

If one could measure all final state momenta in the 3-body decay, 

Y (l{Pq)Q{Pq)x{k), (7) 

where y = ^ or (^(i), and x = i? or -B(i), the spin of Y can be determined by the 2-D 
invariant mass distribution dN^^ecny / df^lgdm'^^ for m^^ = {pg + PqY and m^^ = {pg + fc*''"*^)^ 
or {pq + k^'^'^^)'^. However, as the true WIMP momentum is not available, one could have 



only the m^g-distribution, dNdecay / drri^g = J drn^^dNdecay / drnggdrrig^. In 



, it was found 
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that the SUSY and UED mgg-distributions show a difference in small m^^ limit, however 
this difference might be difficult to be seen in the real data unless the mass ratio my^'^/m^''"'^ 
is quite large, e.g. bigger than 7 or 8. 




m^j (SUSY) [GeV^] fii^j (UED) [GeV^] 



FIG. 2: SUSY and UED Dalitz plots of rn^q and m^^ at parton level for m^^y = m^'^y and very 
large luminosity. 

On the other hand, with the MAOS momenta , we can do a much better job as the 
event distribution dNevent / drn'^^drhg^ is available, where 

= (Pi + ^^)' (P^ + ^^)'- (8) 

In Fig. [21 we depict this MAOS distribution including both and k~ for the SUSY SPS2 
point with m-^ y = m^^'y and its UED equivalent in the ideal limit of no combinatorial 
error and very large luminosity^. The results show a clear difference, with which one can 
distinguish SUSY and UED unambiguously. In fact, these MAOS distributions reproduce 
excellently the shape of the true invariant mass distributions dN^ec&j / drn? drn? constructed 

n 

with the true WIMP momentum k^^^'^ 

To see the feasibility of the MAOS reconstruction in realistic situation, we analyzed the 
event sets of the same SUSY and UED points, but now with the integrated luminosity 
C = 300 fb"^ To suppress the backgrounds, we employed an appropriate selection cut 
commonly taken for new physics events at the LHC, and adopted the hemisphere method 

^ In fact, one can use only or only for the MAOS invariant mass distribution, and still finds the same 
shape of distribution. 
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to deal with the combinatorial problem. We also included the smearing effects on the 
momentum resolution. The results for m^^Y = m^^^ are depicted in Fig. [3l which still 
shows a clear difference between SUSY and UED. 




ml J (SUSY) [GeV^] ^ (UED) [GeV^] 



FIG. 3: SUSY and UED Dalitz plots of m^^ and m^^ for m^,y = m^^^^ and C = 300fb"\ including 
the combinatorial error and smearing effects under a proper event cut. 

A nice feature of the MAOS reconstruction is that it does not require any pre-knowledge 
of the mother particle and WIMP masses. To see if the spin measurement is possible without 
a good knowledge of new particle masses, we have repeated the analysis for a wide range of 
mj(, y. We then found that the basic feature of distribution is retained even when the trial 
masses used for the reconstruction are very different from the true masses, and SUSY and 
UED can be clearly distinguished in all cases [ll]. We depict the result for one example in 
Fig. m which is for the case of = and my = M^2^{m^ = 0), and the result includes 
the combinatorial error and detector smearing effects under a proper event cut. 

Once the MAOS WIMP momenta are available, the MAOS momenta of mother particles 
can be reconstructed also. One can then investigate the production angular distribution of 
mother particle pair in their center of mass frame, which may provide information on the 
spin of mother particle. As an example, we have considered the Drell-Yan pair production 



of Y (= the slepton / or the KK lepton [q], 

qq ^ ^Y + Y^ l^{Vi)x{h) + r(p2)x(A:2), (9) 

and generated the events for the SPSla SUSY point {mf^^ = 143 GeV, m^J^^'^ = 96 GeV) 
and its UED equivalent. 




FIG. 4: SUSY and UED Dalitz plots of m^^ and for = 0, my = M^f'^im^ = 0) and 
C = 300 fb~^, including the combinatorial error and smearing effects under a proper event cut. 

As slepton is a scalar particle, the angular distribution is proportional to 1 — cos^ 6*, where 
6* is the production angle with respect to the proton beam direction. On the other hand, 
the corresponding Drell-Yan production of KK leptons shows the characteristic distribution 
of spin-half particles, which is proportional to 1 + cos^ 6* [Ef^ — mf^ ) / [Ef^ + mf^ ) . One may 
then examine the lepton angular distribution in the center of rapidity frame of //, which would 
reflect the qualitative feature of the above 6'*-dependence of the mother particle distribution 

0. 

Again, we can do a much better job with the MAOS momenta kf {i = 1, 2) as we 
can probe the angular distribution of the mother particle MAOS momenta. To see this, 
we have reconstructed the MAOS momenta, Pi + kf, of the slepton pair and of the KK 
lepton pair, and examined their angular distribution in the center of mass frame while 
including the four different combinations of MAOS momenta, i.e. (kf,k2) with = ±, 
altogether. Since it depends on the longitudinal boost to the center of mass frame, the 
shape of angular distribution is somewhat sensitive to the trial mother particle and WIMP 
masses. To minimize this sensitivity, we have chosen my = M^2^{m^) and imposed the 
event selection cut including only the top 10% of the events near the end-point of Mt2- 
We also included the detector smearing effect on the lepton momentum resolution. In Fig. 
|5^, we depict the resulting SUSY and UED angular distributions for m^^y = m^^y and 
C = 300fb~\ and compare them to the angular distributions obtained from the true WIMP 
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momenta. The result shows that the MAOS angular distribution excellently reproduces the 
true production angular distribution, with which one can clearly distinguish SUSY from 
UED. 

To see how much the distribution shape is sensitive to the trial masses, we repeat the 
analysis, but now with = and my = M^2^{m^ = 0). The result depicted in Fig. 
[5h> shows that the distribution has basically the same shape as the case (Fig. [5^) with 



(energies smeared) 
— SL[SY(tiue) 
• SUSY (MAOS) 
-■- UfiD(triie) 
i UED (MAOS) 


-0^2 0.2 0.4 

ie (in_^ = m™) 



s 

o 

n 160- 



(energies smeared) 
— SUSY (true) 
• SUSY (MAOS) 
-■- UfiD(triie) 
i UED (MAOS) 



-0.8 -0.6 ^-0.4 -0.2 0.2 04 0.6 
COS ( m = Mjj'* {m = 0), m = 0) 



FIG. 5: Slepton and KK-lepton production angular distributions constructed with the true 
momenta and the MAOS momenta for C = 300 fb^^: (a) m^^Y = rn*^^ and (b) = 0, 
my = M^^'^im^ = 0). 

To conclude, we have proposed a scheme, the MT2-assisted on-shell (MAOS) reconstruc- 
tion, to assign 4- momenta to the WIMP pair in generic new physics event of the type ([T]). 
Introducing a trial WIMP mass which is common to the whole event set under consideration, 
the transverse MAOS momenta of each event are determined by the transverse momentum 
that gives the event variable Mt2- On the other hand, the longitudinal MAOS momenta 
are determined by the on-shell condition defined with a trial mother particle mass which 
is again common to the whole event set. With these MAOS WIMP momenta, one can 
construct various kinematic variable distributions which would not be available before, and 
extract information on new particle properties from those distributions. In this paper, we 
considered an application of the MAOS reconstruction to the 3-body decays of the pair- 
produced gluinos and also to the 2-body decays of the Drell-Yan produced charged slepton 



9 



pair to discriminate them from their UED equivalents. Our analysis suggests that the MAOS 
reconstruction of WIMP momenta provides a powerful tool for spin measurement, which can 
be viable even when a good knowledge of the new particle masses is not available. In the 

n 

forthcoming paper [1J.|, we will provide a more extensive discussion of the MAOS recon- 
struction and its applications, including some other processes to measure the new particle 
properties. 
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